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The collision of two intense, low-frequency laser beams is considered. The e _ e + pairs created 
in this field are shown to exhibit recollisions, which take place at high energy accumulated due to 
the wiggling of fermions. The resulting e“e + annihilation produces high energy photons, or heavy 
particles. The coherent nature of the laser field provides strong enhancement of the probability of 
these events. Analytical and numerical results are outlined. 
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A strong static electromagnetic field can create 
electron-positron pairs as was shown by Sauter and 
Schwinger [1, 2]. Similarly, the phenomenon of the pair 
creation takes place in slow-varying fields, as was demon¬ 
strated in [3-7]. This letter elucidates the role of recolli¬ 
sions of the e _ e + pairs created in slowly oscillating fields. 
We show that in a standing wave produced by two collid¬ 
ing laser beams the creation and the following wiggling of 
electrons and positrons in the laser field pave the way for 
their annihilation with high energy gamma-production, 
or production of heavy particles. The probability of these 
phenomena is found to be greatly enhanced by the coher¬ 
ent nature of the laser field. 

The e _ e + production is exponentially suppressed for 
fields below the Schwinger limit £ c = m 2 /e (h = c = 1, if 
not stated otherwise). However, several laser facilities are 
aiming to produce fields within an order of magnitude of 
£ c in the next decade [8-10]. There were also considered 
several ways in which the probability of e _ e + creation 
may be boosted. Refs. [11, 12] discussed manipulation 
of the laser pulse shape and polarization. Other studies 
investigated the dynamical Schwinger mechanism [13], 
in which a weak field is used to lower the barrier for 
the strong Held tunnelling process. Recent work by Di 
Piazza et. al. [14] claims that by tuning the frequency 
of the weak Held, e - e + creation may be observable with 
presently available technology. Meuren et. al. [15, 16] 
considered collision of a single photon with a strong laser 
pulse when the rate of e _ e + creation is enhanced by the 
high photon energy. 

The recollision, or antenna mechanism as originally 
posed in [17], is a three step process. First, an e“e + 
pair is created through the Schwinger mechanism. Sec¬ 
ond, e~ and e + oscillate in the laser field, and third, they 
recollide. These recollisions violate the adiabatic condi¬ 
tion, and hence are exponentially enhanced compared to 
conventional adiabatic processes. In Ref. [17], the recol¬ 
lision of e _ e + pairs created by a heavy nucleus exposed 
to the laser beam were considered, while Meuren et. al. 
[15, 16] considered recollisions in the field of a strong 
laser pulse incident on a single photon. In the problem 
discussed in the present paper, when only the laser fields 


operate, the antenna mechanism reveals a new interesting 
feature, which strongly enhances the collision probability. 
The effect originates from the fact that the velocities of 
the created e~ and e + are correlated along the direction 
of the laser electric field. 

The low energy analogue of the antenna process is the 
atomic antenna mechanism proposed in [18], and devel¬ 
oped quasiclassically by Corkum, and Lewenstein et al. 
[19, 20], in which high-harmonic radiation, multiple ion¬ 
ization and above threshold ionization take place. This 
is currently an area of active research [21, 22], 

Consider the collision of two laser beams represented 
by two plane waves propagating head-on along the in¬ 
direction with the same frequency, intensity, and linear 
polarization along the z-axis. Assume that their phases 
are tuned to produce the standing wave in which the total 
electric field is E = £ cos ut cos kx. We presume that 
the electric field £, frequency w, and adiabatic parameter 
7 = m<jj/(e£) are small, £ < £ c , lo <C m, 7 <C 1; the last 
two inequalities are essential to the e _ e + pair acquiring 
high energies. 

We evaluate the probability W of e“e + creation using 
the Keldysh-type approach [7], in which 

Wa^/|A( P )| 2 ^ . (1) 

Here the integration and summation run over the electron 
momentum p and its spin projection s at the moment of 
creation; one finds that the positron momentum equals 
p while the projections of the electron and positron 
spins along p are opposite. For the amplitude A(p) of 
the pair creation in the relativistic region p 2 m 2 we 
find 


A(p) = A 0 exp(jS'tun(p)) , 

S»(P) = S(- 2 ^i + 7 y l ) 



( 2 ) 

( 3 ) 


where S tAm (p) is the tunneling action and F = e£. Equa¬ 
tions (l)-(3) show that typical electron momenta along 
the electric field greatly exceed perpendicular momenta, 
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p z oc F 1 / 2 /7 p± oc F 1 / 2 . The created fermions there¬ 
fore propagate at small angles to the electric field, which 
complies with expectations based on physical reasons. 
Observe also that (l)-(3) dictate that the pair is cre¬ 
ated mostly in the vicinity of points where the amplitude 
of the electric field is maximal, uit = tv n, ojx = 7 rm, for 
integer n, to. The magnetic field in the standing wave, 
By = £ sinwf svnkx, is small near these points. Cor¬ 
respondingly, on most classical trajectories the fermions 
are exposed only to a small magnetic field and an almost 
homogeneous electric one. Hence below we neglect the 
magnetic field and treat the electric one as homogeneous. 

We aim to calculate the probability of the annihilation 
of the created pair, e _ e + —> 2y, which occurs over a short 
time interval compared to the laser period, u)5t co \ <C 1. 
As a result, the annihilation amplitude can be factorized 
into the matrix element describing the annihilation of 
two fermions with energies e, and the amplitude describ¬ 
ing the fermion propagation to the point of collision. A 
similar factorization takes place in the known problem of 
the positronium annihilation, see e.g. [24]. This analogy 
allows us to write the following expression for the con¬ 
version coefficient, which equals the probability of e - e + 
annihilation, 

T=§/,**(0)| 2 M*)*. ( 4 ) 

Here <j(e) is the cross section of the fermion annihilation 
in the vacuum and v is the electron speed at the moment 
of collision, while <I> 2 £ (ri 2 ) is the wave function of the 
fermion pair. It appears in (4) at zero separation between 
fermions r \2 = 0 to allow the annihilation to take place. 

To evaluate 'f> 2 £ (0) we first construct the time depen¬ 
dent wave function of the fermion pair 'k(ri 2 , t), and then 
examine its Fourier transform. The function 4 '(ri 2 ,t) 
needs to take into account that the pair can be created 
with a given initial electron momentum p; the corre¬ 
sponding amplitude is presented in (2). The propagation 
amplitudes of the fermions are given by their single par¬ 
ticle wave functions ip p \v,t) and where the 

supercript (=p) marks the electron and positron states. 
The pair creation and its following propagation are con¬ 
sequent events. Hence the wave function reads 

$(ri 2 ,t)=P J A(p)V>^" ) ( r i,i)V’^ ) ( r 2 ,i)^ 3 , (5) 

where the integration takes into account interference 
of events with different momenta, and the operator P 
chooses the necessary opposing spin states projections 
of the pair. We find the single-particle wave functions 
using the adiabatic approximation, which is applicable 
since to <C m. The electron wavefunction reads 

^p _) ( r , t) = Xp^t) exp (*( P(t) r ~ J £ { t ') dt ')) > (6) 


where Xp is the usual spinor that describes the elec¬ 
tron propagation with the given momentum p in the 
vacuum. In the time-dependent electric field this mo¬ 
mentum, as well as the energy, become slowly varying 
functions, p (t) = p+(F/w) sinw t, e(t) = (p 2 (t)+m 2 ) 1 ^ 2 . 
The integral in the phase factor, f* e(t')dt', is an impor¬ 
tant and conventional feature of the adiabatic approxi¬ 
mation. The positron wavefunction is constructed similar 
to ( 6 ), with momentum — p(t). 

Substituting Eqs.(2), ( 6 ) into the wavefunction (5) we 
find 

1'(ri2 ,t)=A 0 /epexp(^(p,ri 2 ,t))^, (7) 

S(p, ri 2 , t) = p (t) ■ r i2 - 2 J s(t')dt' + S tu n (p) , ( 8 ) 

where £ P = is the spin part, S(p,r 12 ,t) 

can be considered a classical action, which incorporates 
the tunneling action that describes the pair creation, 
and the action responsible for its subsequent propaga¬ 
tion. The constant Aq here is defined by normalization 
conditions. Since m/w 1, the integral in (7) may be 
evaluated by the saddle point method. The location of 
the saddle point is defined by V p S'(p, ri 2 , t) = 0, in which 
one recognizes the Hamilton-Jackoby condition for a clas¬ 
sical trajectory. This implies that the saddle-point cal¬ 
culations we fulfill amount to a semiclassical description 
of the problem. 

For the case of ri 2 = 0 we are interested in, the rel¬ 
evant classical trajectory has a simple, appealing form. 
Using (3) and ( 8 ) we find that the electron and positron 
velocities have opposite directions and are strictly aligned 
along the electric field, so that the electron and positron 
coordinates can be taken as z(t) and —z(t) respectively. 
The saddle point condition dictates 

z(t) = f v(t') dt' + e(0)/F = 0. (9) 

Jo 

The velocity equals v(t) = p(t)/(p 2 {t) + to 2 ) 1 / 2 , where 
the lower limit of integration in (9) indicates that we 
consider the pair created at ut « 0. We see that at the 
moment of creation the fermions are separated by the 
distance 22 ( 0 ) = 2 e( 0 )/F « 2 |p 2 |/F, which is similar to 
what takes place in the static electric field. The related 
term in the middle expression in (9) was derived from (3) 
neglecting corrections oc y 2 . The last identity in Eq.(9) 
states that the necessary condition ri 2 = 0 takes place at 
the nodes of the function z{t). Since p(t) oscillates, there 
may exist several nodes of z(t), i. e. several collision 
times. From Eq. (9) we find that these nodes appear 
only when p z < 0. Three such nodes are present in Fig. 
1 , where the trajectory z(t) is shown for 7 = 10“ 15 with 
p z = — 0 . 9 F 1 / 2 / 7 . These nodes are enumerated by n = 
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The integral in (11) may be evaluated explicitly, which 
would in turn allow us to calculate the conversion coef¬ 
ficient T in (4) numerically. We elide these details here, 
instead using the wavefunction itself to obtain a reliable 
analytic estimate for T. A simple estimation for T may 
be obtained via classical arguments. One takes into ac¬ 
count that the e + e“ —> 2 7 cross section is of the order 
of a ~ 7rr 2 , see [23], while the typical separation of the 
wiggling fermion pair is R ~ v/to ~ 1/w, because the 
fermion velocity is v ss 1. This leads one to the following 
estimate for the conversion coefficient 


FIG. 1. The electron classical trajectory z(t), 7 = 10 1 ' 8 , 
p z = — 0.9F 1 / 2 /'y. The nodes of z(t) indicate the moments of 
time where e“e + collisions may occur. 


0 , 1, 2 below. 

Evaluating the integral in (7) at the saddle point (9) 
we find 


tf(0,i)=£ p $(0,i)) (10) 

1/2 F 3/4 3?r 

*) ~ E 27/4 Dtp) exp (* (5(p ’ °’- tO 

n> 0 


D{p) 


2n(2?i + u(O)) 1 / 2 In n > 1 

(l + F 1 (3M))( 7 2 y 3 ( 7 ^)+^(0)) 1 / 2 , n = 0 


Tclas ~ oc ~J Uj2f7 ■ ( 12 ) 

However, this result grossly underestimates the probabil¬ 
ity of the recollision because it does not take into account 
that the e _ e + pair propagates predominantly along the 
electric field. We will see that this important property 
stems from the quantum nature of the problem. To jus¬ 
tify this claim we rewrite the conversion coefficient from 
Eq. (4) as follows 

T = 4(w> J l<M0)| 2 g. (13) 

This equality can be considered as the definition of the 
averaged value (va) for the quantity |v|a(e). Using Par- 
seval’s theorem, we can rewrite (13) further 


where Y a (x) = \x — sin (f)\~ a d(f> and n is the number of 
oscillations the fermions undergo before annihilation, see 
Fig. 1 in which the nodes of z(t) correspond to n = 0,1, 2. 

The time t in Eq.(10) is the moment of the collision, 
wlreres p = p z is the electron momentum at the moment 
of creation. It should be chosen in such a way that the 
condition z(t) = 0 is satisfied, which makes p a function 
of t and n. Generally, $(0,/) is a sophisticated function 
of t, which distinguishes the present problem from the 
positronium annihilation where the fermion wavefunction 
exhibits simple oscillations oc exp (—iEp s t), Ep s is the 
total positronium energy. However, the complicated t- 
dependence of the wave function reveals itself only at 
large intervals of time St oc 1/w, while the annihilation 
is a rapid process. Hence, the only contribution relevant 
to the annihilation process originates from the part of 
wave function, which at the moment of collision oscillates 
with the frequency that matches the frequency 2e of the 
photon pair in the final state. To extract this term we 
need to make the Fourier transform 

< f ) 2 e( 0 ) = J exp(— 2iet)$(0,t) dt . (11) 

The spin factor £ p , which is present in (10), is omitted 
here because, following conventional formalism, it is in¬ 
cluded in the definition of the cross section a(s ), which 
appears in (4). 


T = 2 {va) 


|$(0, t)\ 2 dt. 


(14) 


Change now the measure of integration from dt to d\p\. 
The necessary derivative p extracted from the known 
classical trajectory (9), reads |p| = \v\FD 2 (p)/(2n + 
u(0)). Calculating |<I>(0,f)| 2 we presume that the annihi¬ 
lation processes, which take place at different n do not 
interfere with each other since they are separated by large 
intervals of time. To simplify the discussion, we also ne¬ 
glect the term with n = 0 because careful analysis, which 
will be presented elsewhere, shows that numerically its 
contribution is less than that of the n = 1 term. As a 
result, from (10), (14) we derive 


T 


quant 


7 F 1 / 2 {va) ^ 1_ f°° ( 

2 9 / 2 ln 2 (7r/ 7 ) ^ n 2 J Q P V 

^ tt 2 1 F 
192 He 2 ln 2 (7r/ 7 ) 


777V \ d\p\ 
2 F ) M 

(15) 


The limits of integration over \p\ here take into account 
that the collision occurs only for negative values of p. In 
the last step in (15) we presumed that the collision veloc¬ 
ity is large, i>| « 1, which is an adequate approximation. 

Equation (15) presents an appealing, important result. 
The Planck constant there explicitly shows its quantum 
origins, which contrasts the classical nature of Y c i as in 





4 


(12). In the static limit 7 = 0 the annihilation becomes 
impossible, T = 0, because the static electric field can 
but only further separate the created fermions. Equation 
(15) complies with this condition, though for 7 —► 0 the 
term independent of the cross section decreases slowly, 
as oc ln _ 2 ( 7 r/ 7 ). 

Most importantly, our quantum result (15) greatly ex¬ 
ceeds the simple classical estimate ( 12 ) 


- quant 


m 


OC 


Yclas 7 2 ln 2 (7r/7) F 


> 1, 


(16) 


where we presumed that a « (va). Observe a large fac¬ 
tor 7 “, which appears here (the proposed developments 
of laser facilities would make 7 small, 7 ~ 10~ 6 ). The 
arguments presented above show that this impressive en¬ 
hancement stems from the clear physical reason. In the 
coherent laser field the e _ e + pairs are produced with mo¬ 
menta, which are strongly aligned with the electric field. 

To study T quant in more detail we take the two-photon 
annihilation cross section for polarized fermions found 
in [23], which for all energies can be approximated as 
|n|(j(e) s=s 7 r?’ 2 (?n /£) 2 with accuracy > 80%. The typical 
electron energy in the laser field is e oc m/ 7 . Therefore 
(va) ~ nr'l'y 2 , which yields 


tquant 


7 r 3 r 2 7 2 F 
192 ln 2 (7r/7) 


(17) 


This clear and attractive formula is one of the important 
results of this work. We verified its validity by full scale 
numerical calculations based on the formalism developed 
in Eqs.(4), (10), and (11). Our preliminary results are 
shown in Fig. 2. Note that data derived from numerical 
analysis and from Eq.(17) are very close for a wide range 
of 7 , though the distinction grows for larger 7 . 

Figure 2 also shows that at F ss m 2 the maximal value 
of the conversion coefficient is T ss 10 -7 . This means 
that each e~ e + pair created in the laser field can produce 
~ 10 - 7 high energy photons. It should be emphasized 



FIG. 2. The number of photons per e“e + pair T versus 7 = 
muj/F for F = 0.8 m 2 ; dotted line - Eq. (15), solid line - 
numerical results derived from (4), (10), (11). 


that this is a very encouraging result, as this probability 
of the photon production by several orders of magnitude 
exceeds the one predicted for the collision of a laser beam 
with the nucleus [17]. As was mentioned, the reason for 
the found enhancement stems from the coherent nature 
of the laser field. 

We can describe similarly the production of heavy, high 
energy particles, for example muons, triggered by the 
annihilation of e _ e + pairs created in the standing laser 
wave. To make this phenomenon possible 7 should be 
sufficiently small to ensure that the wiggling energy of 
the e~ and e + to exceed the mass M of a heavy particle, 
e ~ m /7 > M. For 7 ~ 10 ~ 6 this wiggling energy is 
quite large, ~ 1 TeV, as one derives from Fig. 3, where 
|$ 2 £ ( 0)| 2 measures the probability that the e _ e + collision 
takes place with the energy e. The antenna mechanism 
considered gives a huge, exponential enhancement for 
such processes compared with the direct, Schwinger-type 
production of heavy particles. The probability for the lat¬ 
ter is oc exp(— ttM 2 /F), whereas the antenna mechanism 
gives T oc exp(—irm 2 /F), where m is the electron mass. 
According to Eq.(15) the main suppression for T orig¬ 
inates from the cross section, which decreases with en¬ 
ergy. However, it is a power-type decrease, while the gain 
in probability is exponential, oc exp ( 7 t(M 2 — m 2 )/F). 
Similarly, by comparing the production of high energy 



2s (TeV) 

FIG. 3. Fourier component of the wave function, |4>2e(0)| 2 
(arbitrary units), numerically calculated from Eq.(ll) for n = 
0 versus the electron energy e for 7 = 2 - 10 ~ 6 ; bottom to top: 
F/m 2 = 0.0125, 0.025, 0.05, 0.1, 0.2, 0.4, 0.8. 

photons by an independent fermions e^, which radiate 
due to their oscillation in the laser field, and the pho¬ 
ton production via e _ e + annihilation, considered in this 
work, we find that later process is exponentially enhanced 
compared to the former one. 

Summarizing, it is demonstrated that in the foresee¬ 
able future high energy photons and massive particles 
with energies above ~ 1 TeV can be produced using the 
intense colliding laser beams. 

This work was supported by the Australian Research 
Council. 
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